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I. INTRODUCTION 

The theory that describes the fundamental par- 
ticle interactions is called the standard model 
(SM). The standard model is a gauge field the- 
ory that comprises the Glashow-Weinberg-Salam 
(GWS) model Jl|-[| of the weak and electromag- 
netic interactions and quantum chromodynamics 
(QCD) the theory of the strong interactions. 

The discovery of the W @f| and Z ||[lO| bosons 
in 1983 by the UA1 and UA2 collaborations at the 
CERN pp collider provided a direct confirmation of 
the unification of the weak and electromagnetic in- 
teractions. Experiments have been refining the mea- 
surements of the characteristics of the W and Z 
bosons. The total decay width of W boson, Tw, is 
given in the SM in terms of the masses of the gauge 
bosons and their couplings to their decay products. 

In pp collisions, W bosons are produced by pro- 
cesses of the type ud or ud — > W, followed by sub- 
sequent leptonic or hadronic decay: W — > £v or 
W — > q'q, where i = e,(x,T, and q' or q represent 
one of the quarks u, d, c, s or b (but not t since top 
quark is heavier than the W boson). 

At lowest order in perturbation theory, the SM 
predicts the partial decay width T(W — ► ev) of W — ► 
ev to be T(W -> ev) = g 2 M w /48tt 0. Including 
radiative corrections, this can be rewritten as: 
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where Gf/V% — g 2 /8M w , g is the charged current 
coupling, and Mw is the mass of the W boson. The 
SM radiative correction, Ssm, is calculated Q to be 
less than By using the experimental values of 
Gf (measured from muon decay |l3|]) and Myy (mea- 
sured at the Fermilab Tevatron collider [jl4],[l5) and 
LEP2 JT§-0), the predicted partial width is Q 
T(W -> ev) = 226.5 ± 0.3 MeV 

A W boson has three leptonic decay channels 
and two dominant hadronic decay channels, W — > 
ev,/j,v,TV and qq', where q is u or c, and q' is 
the appropriate CKM mixture of d and s. Other 
hadronic decay channels arc greatly suppressed by 
CKM off-diagonal matrix elements. Considering 
the three color charges for quarks, these nine lep- 
tonic and hadronic channels yield a total width of 
« 9r(14 7 — > ev). Including QCD corrections, the 
leptonic decay branching ratio is B(W — ► ev) — 
l/{3 + 6[1 + a s (M w )/n + 0(a 2 )}}, leading to the 
SM prediction for the full width of the W boson Q 
of T w = 2.0921 ± 0.0025 GeV. 

Historically, the accurate determination of the 
width of the W boson was available through an indi- 
rect measurement using the ratio 1Z of the W — > ev 
and Z — > ee cross sections: 



a(pp — > Z 4 
a w Br(W 
~a~z~ ' Br(Z 



X) ■ Br(Z 
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A measurement of TZ, together with a calculation ECj 
of the ratio of production cross sections aw/cz and 
the measurement of the branching faction Br(Z — > 
ee) = T{Z ee)/T(Z) from the CERN e+e" col- 
lider (LEP) [^l|, can be used to extract the W boson 
leptonic branching ratio Br(W — > ev) — T(W — > 
ev)/T(W), which, in turn, yields the full width of 
the W boson from calculated partial decay width 
T(W — > ev). Thus, in this indirect measurement, 
calculations of aw /az and the partial width T(W — > 
ev) yield Tw in the context of the SM. This method 
was first used by the UAl H and UA2 f|| collab- 
orations. More recently, the CDF (24) and D0 Q 
collaborations obtained T w = 2.064±0.084 GeV and 
Tw = 2.169 ± 0.079 GeV, respectively, using this 
technique. 

The value of Tw can also be obtained from the 
line shape of the transverse mass tot of the W bo- 
son, because the Breit-Wigner (width) component of 
the line shape falls off more slowly at high tht than 
the resolution component does |p^| . The transverse 
mass is given by 



m T = J2E^E^[1 - cos( 



(3) 



where Ej. and Elf are the transverse energies, and 
<f> e and <ff are the azimuthal angles, of the elec- 
tron and neutrino, respectively. The transverse mass 
has a kinematic upper limit at the value of Mw, 
and the shape of the tot distribution at this up- 
per limit, called the "Jacobian edge," is sensitive 
to T w @- Using this technique, the CDF col- 
laboration reported J27J a measurement of Tw = 
2.05 ± 0.10(stat.)±0.08(syst.) GeV. Figure | shows 
the rriT spectrum shape expected for different val- 
ues of Tw and indicates the sensitivity of the tail of 
the transverse mass distribution to Tw- Clearly, the 
effect is greatest in the region above mw- 

The direct measurement of Tw complements the 
indirect measurement through 1Z in several ways: 

• Theoretical inputs for aw/az and T(W — > 
ev), which may be sensitive to non-SM cou- 
pling of the W boson, are not needed. 

• The direct measurement explores the region 
above the W boson mass pole, where possible 
new phenomena such as an additional heavy 
vector boson (W) can contribute. 

• It is desirable to have more than one method 
of measuring a given property. The sources 
of systematic errors in the two methods are 
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FIG. 1. Monte Carlo simulations of the transverse 
mass spectrum for different W boson widths. The se- 
lections, E T (e) > 25 GeV and E T {v) > 25 GeV, are 
applied to MC sample. The circles show the spectrum 
for IV = 1.60 GeV, the squares for IV = 2.10 GeV, 
and triangles for Tw = 2.60 GeV. Distribution are nor- 
malized arbitrarily in the transverse mass region shown. 



different, and the direct method will be im- 
portant when the measurement through 1Z be- 
comes limited by systematic uncertainty. 

The paper is organized as follows. In Sec. II, we 
give a brief description of the D0 detector. Parti- 
cle identification and event selection are discussed 
in Sec. III. The analysis procedure, including back- 
ground estimation and Monte Carlo simulation, is 
described in Sec. IV, and the conclusions are pre- 
sented in Sec. V. For more detailed information on 
this analysis, see Ref. 



II. THE D0 DETECTOR 



A. Experimental Apparatus 



The D0 detector 1 30 1 comprises three major sys- 
tems. The innermost of these is a non-magnetic 
tracker used in the reconstruction of charged parti- 
cle tracks. The tracker is surrounded by central and 
forward uranium/liquid- argon sampling calorime- 
ters. These calorimeters are used to identify elec- 
trons, photons, and hadronic jets, and to reconstruct 
their energies. The calorimeters are surrounded by 
a muon spectrometer used in the identification of 
muons and the reconstruction of their momenta. We 



use a coordinate system (p, 9, <fi) where p is the per- 
pendicular distance from the beam line, 9 is the po- 
lar angle measured relative to the proton beam di- 
rection z, and <f> is the azimuthal angle. The pseudo- 
rapidity rj is defined as — ln(tan |). For this analysis, 
the relevant components are the tracking system and 
the calorimeters. 

The central tracking system provides a measure- 
ment of the energy loss due to ionization (dE/dx) 
for tracks within its tracking volume. This infor- 
mation is used to help distinguish prompt electrons 
from e + e~ pairs due to photon conversions. 

The structure of the calorimeter has been op- 
timized to distinguish electrons and photons from 
hadrons and to measure their energies. It is 
composed of three sections: the central calorime- 
ter (CC), and two end calorimeters (EC). The rj- 
coverage for electrons used in this analysis is | rj |< 
1.1 |^9| in the CC region, which consists of 32 <fi mod- 
ules. The calorimeter is segmented longitudinally 
into three sections, the electromagnetic calorimeter 
(EM), the fine hadronic calorimeter (FH), and the 
coarse hadronic calorimeter (CH). The EM calorime- 
ter is subdivided longitudinally into four layers 
(EM1-EM4). The first, second and fourth layers 
of the EM calorimeter are transversely divided into 
cells of size Ar/xA0 = 0.1x0.1. The electromagnetic 
shower maximum occurs in the third layer, which is 
divided into finer units of 0.05 x 0.05 to improve the 
measurement of the shower shape and spatial resolu- 
tion. There are 16 FH modules and 16 CH modules 
in <f>. The fine hadronic calorimeter is subdivided 
longitudinally into three fine hadronic layers (FH1- 
FH3), and there is only one coarse hadronic layer. 



B. Trigger 

The D0 trigger has three levels, each applying in- 
creasingly more sophisticated selection criteria to an 
event. The lowest level trigger, Level 0, uses scin- 
tillation counters located on the inner faces of the 
forward calorimeters to signal the presence of an in- 
elastic pp collision. Data from the Level counters, 
the calorimeter, and the muon chambers are sent to 
the Level 1 trigger, which provides a trigger on total 
transverse energy {E^), missing transverse energy 
(fx)-: of individual calorimeter towers, and/or 
the presence of a muon. These triggers operate in 
less than 3.5 fj,s, the time between bunch crossings. 
Some calorimeter and muon-based triggers require 
additional time, which is provided by a Level 1.5 
trigger system. 

Level 1 (and 1.5) triggers initiate a Level 2 trig- 
ger system that consists of a farm of microproces- 
sors. These microprocessors run simplified versions 
of the off-line event reconstruction algorithms to se- 
lect events of interest. 
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III. PARTICLE IDENTIFICATION AND 
EVENT SELECTION 

This analysis relies on the D0 detector's abil- 
ity to identify electrons and to associate the unde- 
tected energy with neutrinos. We use both W — » ev 
and Z — > e + e~ candidate samples for this analy- 
sis. The W boson candidate sample provides the 
signal events, while the Z — > e + e~ candidate sam- 
ple is used to calibrate both the data and the Monte 
Carlo (MC) simulation. Candidate W and Z events 
are identified by the presence of an electron and a 
neutrino, or by the presence of two electrons with 
an invariant mass consistent with the mass of the 
Z boson, respectively. Electrons from W and Z bo- 
son decays typically have large transverse energy and 
are isolated from other particles. They are associ- 
ated with a track in the tracking system and with a 
large deposit of energy in one of the EM calorime- 
ters. Neutrinos do not interact in the detector, and 
thus create an apparent transverse energy imbalance 
in an event. For each W boson candidate event, we 
measure the energy imbalance in the plane trans- 
verse to the beam direction ($ T ), and attribute this 
to the neutrino. The following sections provide a 
brief summary of the procedure pq] used in this 
analysis. 



A. Electron Identification 

Identification of electrons starts at the trigger 
level with the selection of clusters of electromag- 
netic energy. At Level 1, the trigger searches for 
EM calorimeter towers (A</> x A-q = 0.2 x 0.2) with 
signals that exceed predefined thresholds. W bo- 
son triggers require that the energy deposited in a 
single EM calorimeter tower exceed 10 GeV. Those 
events that satisfy the Level 1 trigger are processed 
by the Level 2 filter. The trigger towers are com- 
bined with the energy in the surrounding calorimeter 
cells within a window of x A?/ = 0.6 x 0.6. 

Events are selected at Level 2 if the transverse 
energy in this window exceeds 20 GeV. In addition 
to the Et requirement, the longitudinal and trans- 
verse shower shapes are required to match those ex- 
pected for electromagnetic showers. The longitudi- 
nal shower shape is described by the fraction of the 
energy deposited in each of the four EM layers of the 
calorimeter. The transverse shower shape is charac- 
terized by energy deposition patterns in the third 
EM layer. The difference between the energies in 
concentric regions covering 0.25x0.25 and 0.15x0.15 
in Ar] x A(j> must be consistent with that expected 
for an electron j30j . 

In addition, the electron candidates are required 
to deposit at least 90% of their energy in the EM 
section of the calorimeter and to be isolated from 



other calorimetric energy deposits. To be considered 
isolated, electrons must satisfy the isolation require- 
ment, /i SO < 0.15, where /i SO is defined as: 



fftotal(0.4)-ff E M(0.2) 

£em(0.2) 



(4) 



in which E ioia \(0A) is the total energy, and £?em(0.2) 
the electromagnetic energy, in cones of radius R — 
^(Ar/) 2 + (A0) 2 = 0.4 and 0.2, respectively. This 
enhances the signal expected from isolated electrons 
in W and Z boson decay. 

Having selected events with isolated electromag- 
netic showers at the trigger level, we first define 
"loose" electron for the purpose to study the back- 
ground. Those EM clusters are require to locate 
within the center 80% of a calorimeter module, have 
an associated track in the central tracking volume 
and I r] I < 1.1. To avoid areas of reduced response 
between neighboring calorimeter modules, the az- 
imuthal angle of electrons is required to be at least 
A<fi = 0.10 x 27r/32 radians away from the position 
of a module boundary. We further impose a set of 
off-line tighter criteria to identify electrons, thereby 
reducing the background from QCD multijet events. 
The first step in identifying an electron is to form 
a cluster around the trigger tower using a nearest 
neighbor algorithm. As at the trigger level, the clus- 
ter is required to be isolated (/; so < 0.15). To in- 
crease the likelihood that the cluster is due to an 
electron and not a photon, a charged track from the 
central tracking system is required to point to the 
center of the EM cluster. We extrapolate the track 
to the third EM layer of the calorimeter and cal- 
culate the distance between the extrapolated track 
and the cluster centroid along the azimuthal direc- 
tion (pAcjj) and in the z-direction (Az). The posi- 
tion of cluster centroid is defined at the radius of 
the third EM layer of the calorimeter. The z posi- 
tion of the event vertex is then defined by the line 
connecting the centroid position of the EM cluster 
to associated one in in the central tracking system 
and extrapolated to the beam line. The electron Et 
is calculated using this vertex definition |25|] . The 
variable 



'trk 



V °>0 



) 



Az 



(5) 



where a p( p and a z are the respective track resolu- 
tions, quantifies the quality of the match. A re- 
quirement of <7trk < 5 is imposed on the data. These 
clusters are then subjected to a 4- variable likelihood 
test ]3l| , |32| ]. The four variables are: 

• A \ 2 comparison of the shower shape with the 
expected shape of an electromagnetic shower, 
computed using a 41-variable covariance ma- 
trix [p3[ for the energy depositions in the cells 
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of the electromagnetic calorimeter and the lo- 
cation of event vertex. 

• The electromagnetic energy fraction, defined 
as the ratio of shower energy in the EM section 
of the calorimeter relative to the sum of EM 
energy plus the energy in the first hadronic 
section of the calorimeter. 

• A comparison of the track position to the po- 
sition of cluster centroid, as defined in Eq. |5j. 

• The ionization, dE/dx, along the track. This 
is used to reduce contamination due to e + e~ 
pairs from photon conversions, mainly from 
jets fragmenting into neutral pions. The e + e~ 
pair from photon conversion has a double value 
of dE I dx for a geniue electron due to two over- 
lapping tracks, 

To good approximation, these four variables are in- 
dependent of each other for electron showers. Elec- 
trons that satisfy all above criteria are called "tight" 
electrons. 

Electron energies are corrected for the underlying 
event energy that enter into the electron windows. 
The electromagnetic energy scale is determined in 
the test beam data, and adjusted to make the peak of 
the Z — > e + e~ invariant mass agree with the known 
mass of the Z boson [flj. We found it to be 0.9545 ± 
0.0008. The electron energy scale is discussed in 
detail in Ref. ||. 

B. Missing Transverse Energy 

The primary sources of missing energy in an 
event include the neutrinos that pass through the 
calorimeter undetected and the calorimeter resolu- 
tion. The energy imbalance is measured only in the 
transverse plane because of the lost particles emitted 
at small angles (within the beam pipes). The miss- 
ing transverse energy is calculated by taking the neg- 
ative of the vector sum of the transverse energy in 
all of the calorimeter cells. This gives both the mag- 
nitude and direction of $ T , allowing the calculation 
of the transverse mass of the W boson candidates. 



C. Event Selection 

The W boson data sample used in this analysis 
was collected during the 1994-1995 run of the Fer- 
milab Tevatron collider, and corresponds to an in- 
tegrated luminosity of 85.0 ± 3.6 pb _1 . Events are 
selected by requiring one tight electron in the central 
calorimeter (| rj |< 1.1) J29| with E T > 25 GeV. In 
addition, events are required to have $ T > 25 GeV 
and W transverse momentum px{W) < 15 GeV, 



which is combined transverse momentum of electron 
and $ T (neutrino). After applying all of the de- 
scribed selections, a total of 24487 W boson can- 
didates is selected. There are 24479 candidates in 
the region - 200 GeV, while 8(2) candidates have 
m,T >200(250) GeV. Figure || shows the transverse 
mass distribution of the W — > ev candidates. 

Candidates for the process Z — > e + e~ are required 
to have two tight electrons, each with Et > 25 GeV 
in the CC. The invariant mass of the dielectron pair 
is required to satisfy 60 GeV < m ee < 120 GeV. A 
total of 1997 Z boson candidates is selected. Fig- 
ure H shows the invariant mass distribution of the 
Z — > e + e~ candidates. 




40 60 80 100 120 140 160 180 200 
Transverse mass (GeV) 

FIG. 2. Transverse mass distribution of W — ► ev 
event candidates. 



IV. ANALYSIS PROCEDURE 

In this section, we describe the Monte Carlo sim- 
ulation program used to model the transverse mass 
spectrum. The background from the dominant pro- 
cesses that can mimic the W — > ev signal is also 
estimated. We compare the data with the expecta- 
tion from the Monte Carlo simulation and extract 
the decay width of the W boson using log-likelihood 
fits to the W boson transverse mass distribution. 



A. Monte Carlo Simulation 

The transverse mass spectrum for the W boson 
is modeled in three steps: W boson production, W 
boson decay, and a parameterized detector Simula- 

tion mnii. 
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FIG. 3. Invariant mass distribution of Z — + e e 
events compared to Monte Carlo simulation. The histro- 
gram is the MC and the black dot with error bar is the 
data. The Z — > e + e~ candidate require both electrons 
be in the CC. 



We first simulate the production of the W boson 
by generating its four momentum and other event 
characteristics, such as the z-position of the inter- 
action vertex and the run luminosity. The luminos- 
ity is used to parameterize luminosity-dependent ef- 
fects. To lowest-order, the mass of the W boson 
follows the Breit-Wigner distribution: 



a(Q) = C qq {Q) 



Q 2 



(Q 2 



y + Q^ w /M* w 



(6) 



qq(Q) is called the parton luminosity. To 
qqve ii rr - ev events using 



where Q is the invariant mass of W boson, Myy is 
the pole mass and Yy/ the decay width of the W bo- 
son, and L 

evaluate C qq -(Q), we generate W 
the leading-order resbos Q event generator and 
the different PDF models described in Refs. fl37| , [38| ] , 
and then select the events using the same kinematic 
and fiducial constrains as for the W and Z boson 
data samples. The resulting event distribution is 
proportional to the parton luminosity, which we pa- 
rameterize with the function |39l: 



£qq(Q) = 



(7) 



where (3 is obtained from a fit of the MC events to 
Eq. §. 



The decay of the W boson is simulated in the MC 
and used to calculate the transverse momentum of 
the electron and other decay products. Any radi- 
ation from the decay electron or from the W bo- 
son can bias the measurement and has to be taken 
into account. W — > rv — > evDD events are indis- 
tinguishable from W — > ev and are also included 
in the model, using a branching ratio of Br(r — > 
ew)/[l + Br(r -> evv)} = 0.151. 

Finally, we apply a parameterized detector simu- 
lation to the momenta of all decay products to simu- 
late any observed recoil jets and electron momenta. 
The parameters giving the electron and recoil sys- 
tem response of the detector are fixed using data, 
which include Z bosons and their recoil jets, to study 
calorimeter response and resolution. The response 
to jets and electrons is parameterized as a function 
of energy and angle. Also included in the detector 
parameterization are effects due to the longitudinal 
spread of the interaction vertex and the luminosity- 
dependent response of the detector caused by mul- 
tiple collisions. 

Uncertainties in the input parameters to the MC 
will eventually limit the accuracy of the width mea- 
surement of the W boson. To study the uncertain- 
ties, we allow these input parameters to vary by one 
standard deviation and re-generate the correspond- 
ing transverse mass spectrum. We then fit it with a 
nominal MC template. If the positive and negative 
variations of the width of the W boson with respect 
to a parameter are not symmetric, the larger value 
is used for the uncertainty. This estimation is used 
to estimate the impact of the electron energy reso- 
lution, hadronic energy resolution, electron energy 
scale, hadronic energy scale, dependence on the W 
boson mass, electron angular calibration, and radia- 
tive corrections. Detailed studies of these parame- 
ters can be found in Ref. Jl5[ |. The uncertainties on 
Y\y from the electron energy resolution and scale are 
27 MeV and 41 MeV, respectively. The uncertainties 
from the hadronic energy resolution and scale lead 
to variations in Tw of 55 MeV and 22 MeV, respec- 
tively. The error on the W boson mass of 37 MeV 
has an effect of 15 MeV on Fw The uncertainties 
from radiative decay and electron angular calibra- 
tion correspond to 10 MeV and 9 MeV, respectively. 

Uncertainties on Tw also arise from uncertain- 
ties in the production model and the parton dis- 
tribution functions (PDFs). The uncertainty from 
the former is determined from the upper and lower 
limits p7| of the most uncertain parameter in the 
model. This leads to an uncertainty of 28 MeV 
due to parton luminosity and 12 MeV due to un- 
certainty in the transverse momentum of the W bo- 
son in the model. There are several PDF models 
currently in use. The uncertainty due to variation 
in PDFs is determined by using different PDFs, in- 
cluding MRSA @, CTEQ4M and CTEQ5M 
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and finding the largest excursion from the value of 
IV determined using the MRST PDF set [||, lead- 
ing to a variation of 27 MeV. The value quoted for 
IV is determined using the MRST PDFs. We chose 
MRST so that the results can be consistent with D0 
mass analysis |jl5| . 

B. Backgrounds 

Backgrounds to W — > ev can affect the shape of 
the rriT spectrum and skew the measurement of r^. 
We account for this by estimating the background as 
a function of vtlt and adding this to the tot dis- 
tribution of the W boson from the Monte Carlo. 
The three dominant background sources are multi- 
jet events, Z — > ee, and W — » tv decay products. 
The following describes how the backgrounds are es- 
timated ||§. 

A large potential source of background is due to 
multijet events in which one jet is misidcntified as 
an electron and the energy in the event is mismea- 
sured, thereby yielding large $ T . This background 
is estimated using jet events from data, following the 
procedure called the "matrix method," described in 
Rcf. |25|,|28|,|32|] . The method uses two sets of data, 
each containing both signal and background. The 
first data set corresponds to the W data sample in 
this analysis. The second set contains a different 
mix of signal and background which is obtained with 
loose electron criteria (described in Sec. Ill A). We 
summarize below the essence of this method used to 
estimate the multijet background. 

The number of multijet background (N^ G ) events 
in the tight electron W data sample is given by 



where Ni and N t are the number of events in the 
W boson samples satisfying loose and tight elec- 
tron criteria, respectively. The tight electron ef- 
ficiency, e s , is the fraction of loose electrons that 
pass tight electron criteria, as determined by the 
Z boson sample, where one electron is required to 
pass the tight selection criteria and the other serves 
as an unbiased probe for determining relative effi- 
ciencies. The electron efficiency is obtained to be 
e s = (86.3 ± 1.2)%. The jet efficiency £j is the frac- 
tion of loose "electrons" found in multijet events 
that also pass tight electron criteria. This sample 
is required to have $! T < 15 GeV to minimize the 
number of W bosons contained in it. The results is 
€j = (5.83 ±0.25)%. Once e s and Cj are determined, 
we can extract the background-event distribution. 
The "electron" and "neutrino" transverse momenta 
and energies are used to form the transverse mass, 
and this distribution is shown in Figure ^. The to- 
tal multijet background is estimated to be 368 ± 32 



events in the region m T < 200 GeV, with 25.4 ± 2.2 
events in the range 90 GeV < tut < 200 GeV. 

The background sample is smoothed in the region 
85 GeV < m T < 200 GeV. We fit the distribu- 
tion to an exponential function of the form Jbg = 
exp(ao + a\x + a^x 1 + 03a; 3 ). The fitting parameters 
an, ffli, a>2 and a 3 p3| are used to generate the back- 
ground distribution for the fit to the signal. For bins 
outside the fitted region, we use the original data it- 
self, as shown in Figure [|. 

Another source of background is due to Z — > ee 
events in which one electron is undetected. This 
results in a momentum imbalance, with the event 
now being topologically indistinguishable from W — > 
ev events. This background is also estimated using 
Monte Carlo events. The number of such Z boson 
events present in the W boson sample is calculated 
by applying the W boson selection criteria to MC 
Z — ► ee events generated using herwig Q] and 
processed through a GEANT J45|] based simulation 
of the D0 detector, and then overlaid with events 
from random pp crossings. This is done to simulate 
the effect of the luminosity on the underlying event. 
Out of a total of 8870 Z — > ee events, 48 pass the W 
boson event selection. Normalizing the Monte Carlo 
sample to the size of the data sample for equivalent 
luminosity, we estimate that there are 102 Z — ► ee 
events in the data sample. 



> 10 2 - 
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FIG. 4. The transverse mass distribution for the mul- 
tijet background. The line represents the results of the 
fit described in the text. 

W ' — > tv events in which the r decays into an elec- 
tron and two neutrinos are indistinguishable from 
W — > ev events on an event-by-event basis. Because 
t undergoes a three-body decay, leading to a softer 
electron relative to W — > ev events, the acceptance 
is reduced greatly by the standard Et selection cri- 
teria. The size of this background is small, and it 
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tends to add events with low values of mp. This 
background is determined using the W — > ev Monte 
Carlo, modified to include the decay of the r lepton. 
The events are then passed through the same detec- 
tor simulation used to model the W — ► ev signal. 

The shape and total amount of background af- 
fect the fit used to determine the width of W boson. 
To estimate the uncertainty in Y\y due to the un- 
certainty in absolute background, we scale up (and 
down) the fitted number of background events by an 
amount that corresponds to the total uncertainty 
in the background. This gives an uncertainty of 
15 McV for Tw extracted from the region 90 GeV 
< rriT < 200 GeV. To estimate the uncertainty in 
F\y from the uncertainty in the shape of the back- 
ground spectrum, we perform an ensemble study in 
which background is generated using a multinomial 
distribution. The multinomial distribution is defined 
by: 

eft Ni 

P(N 1 ,N 2> ...,N ch ) = N total l[^. (9) 

i— 1 l ' 

where A^totai is the total number of background 
events, ch is the number of the bins, pi is the origi- 
nal distribution, and Ni is numbers of events in i-th 
bin. The total background iVtotal is kept at its cen- 
tral value, while the number of background events 
in each bin is allowed to fluctuate. The W boson 
width is then recalculated with the new background 
distribution. The variation in Tyy is taken as the 
uncertainty. We found that this is 39 MeV for the 
fitted region of rriT- 



C. Likelihood Fitting 

We generate a set of Monte Carlo tut templates 
with Tw varying from 1.55 GeV to 2.75 GeV at in- 
tervals of 50 MeV. These templates are normalized 
to the number of events in the region of < 200 
GeV. The background distributions of multijet and 
Z — > ee events are added to the templates and a 
binned likelihood is calculated for data. The tut 
bin size is 5 GeV. The fitting region is chosen to 
be 90 GeV < mx < 200 GeV to minimize the sys- 
tematic uncertainty. From the dependence of the 
likelihood on Tw> we obtain the W boson width 
and its error as: T w = 2.23±£^(stat.) GeV. 
The combined uncertainty, taking the statistical and 
systematic uncertainties contribution in quadrature, 
yields the result T w = 2.23±^(stat.) ± 0.10(sys.) 
GeV= 2.23±£$ GeV. The x 2 for the best fit is an ac- 
ceptable 25.9 for 22 degrees of freedom, correspond- 
ing to a probability of 26%. A comparison of the ob- 
served spectrum to the probability density function 
in the fitting region through a Kolmogorov-Smirnov 



test, which compares the observed cumulative distri- 
bution function for a variable with a specified the- 
oretical distribution, yields n — 0.434, which is evi- 
dence of a good fit. 

Figure || shows a fit to the likelihood, which cor- 
responds to a fourth-order polynomial fit that de- 
termines the peak position. Figure ^| shows the mj- 
spectrum for the data, the normalized MC sample, 
and the background. 
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FIG. 5. Results of the log-likelihood fit of the data 
to Monte Carlo templates for different Tw- 

As a consistency check of the fitting method, we 
also determine the W boson width from the ratio of 
the number of events in the fitting region of 90 GeV 
< hit < 200 GeV to the number of events in the 
entire spectrum. This yields = 2.22±0.14(stat.) 
GeV, compared to T w = 2.23±j^(stat.) GeV for 
the independent maximum likelihood fit in the same 
region. All results show good agreement. 

Sources of systematic uncertainties in the deter- 
mination of the W boson width are those that can 
affect the shape of the transverse mass distribution. 
These include the uncertainties from input parame- 
ters to the MC program and from background esti- 
mation. Details can be found in Ref. |^8| . Table [j| 
lists all the sources of systematic uncertainty for the 
decay width of the W boson. 

Comparing to the SM prediction of r(VF) = 
2.0921 ± 0.0025 GeV, we find the difference between 
SM prediction and our meaasurement to be 0.241^ ^ 
GeV, which is the width for the W boson to decay 
into final states other than the two lightest quark 
doublets and the three lepton doublets. We set a 
95% confidence level upper limit on the W boson 
width to non-SM final states. Assuming the uncer- 
tainty is Gaussian, we set a 95% confidence level 
upper limit on the invisible partial width of the W 
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FIG. 6. Comparison of data to the Monte Carlo tem- 
plates for the best fit. The black circles with error bars 
are the data. The solid line of the histogram corresponds 
to the MC templates with V{W) = 2.23 GeV normalized 
to the expected number of W boson events. The shad- 
owed area is the background. 



Source 


Sr w (MeV) 


Hadronic energy resolution 


55 


EM energy scale 


41 


Background ensemble studies 


39 


Luminosity slope dependence 


28 


EM energy resolution 


27 


PDF 


27 


Hadronic energy scale 


22 


Background normalization 


15 


W boson mass 


15 


Production model 


12 


Radiative correction 


10 


Selection bias 


10 


Angular calibration of e trajectory 


9 


Total systematic uncertainty 


99 


Total statistical uncertainty 


+145 




-138 


Total uncertainty 


+176 




-170 



TABLE I. Systematic uncertainties and the total un- 
certainty on the W boson width measurement. 



boson to be 0.59 GeV. Under the assumption that 
there is no correlation between indirect measurement 
and direct measurement of the W boson decay width 
and within the framework of SM, we can combine 
both analysis and obtain Tw = 2.162 ± 0.062 GeV. 
The 95% confidence level upper limit on the invisible 
partial width of the W boson is 0.191 GeV. 



V. CONCLUSIONS 

We have directly measured the decay width of the 
W boson by fitting the transverse mass in W — > eu 
events in pp collisions at 1.8 TeV, and obtain: 

r w = 2.23l°;^(stat.) ± 0.10(syst.) GeV (10) 
= 2.23±g;if GeV. (11) 

This result is consistent with the prediction of the 
standard model. 
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